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Abstract 



This work presents a 2.5-dimensional simulation study of the instabihty of 
current-sheets located in a medium with a strong density variation along the 
current layer. The initial force-free configuration is observed to undergo a two- 
stage evolution consisting of an abrupt regime transition from a slow to a fast 
reconnection process leading the system to a final chaotic configuration. Yet, 
the onset of the fast phase is not determined by the presence of any anomalous 
enhancement in plasma's local resistivity, but rather is the result of a new mech- 



anism discovered in Lapenta (2008) and captured only when sufficient resolution 
is present. Hence, the effects of the global resistivity, the global viscosity and 
the plasma beta on the overall dynamics are considered. This mechanism allow- 
ing the transition from slow to fast reconnection provides a simple but effective 
model of several processes taking place within the solar atmosphere from the high 
chromosphere up to the low corona. In fact, the understanding of a spontaneous 
transition to a self-feeding fast reconnection regime as well as its macroscopic 
evolution is the first and fundamental step to produce realistic models of all 
those phenomena requiring fast (and high power) triggering events. 



Subject headings: MHD - Sun: activity - Sun: chromosphere - Sun: magnetic field 
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Introduction 



Magnetic reconnection is a fundamental mechanism playing a major role both in the 
dynamics of astrophysical objects as well as in the evolution of several structures in the 
solar environment, in planetary magnetospheres and in the heliosphere. Here, indeed, 
a broad range of phenomena undergoes a multi-stage evolution determined by a rapid 
transition from a slow to a fast, and often bursty, dynamics which presents typical magnetic 
reconnection marks as the change in the magnetic field-line connectivity, the formation of 



high- velocity plasma jets, plasma heating, and particle accelerations (Biskamp 1993). 



Magnetic reconnection is believed to be the key process underlying the evolution of 
a number of processes that take place in the solar environment. In particular, it may 
be the responsible for the release of magnetic flux and plasma in the corona and in the 
interstellar medium during solar explosive phenomena. In general, a catastrophic loss 
of mechanical equilibrium of the underlying coronal magnetic structure is required and 
magnetic reconnection plays a role either in triggering the initiating eruption or in the 
aftermath of the event by characterizing the subsequent evolution and related structures 



like post CME's current-sheets (Forbes & Lin 2000 Lin et al. 2004). In the case of a 



flare evolution, the free energy of the solar magnetic field (exceeding the energy of the 
solar atmosphere potential field) is released. This excess energy is connected with currents 
flowing in the corona from the photosphere throughout the chromosphere and the transition 
region, the flare process being the process of rapid variation of these currents. 



From Yohkoh X-ray telescope observations ( Shibata||1996 ), we have evidences consistent 



with the reconnection model of solar flares (Tsuneta 1996). In particular, the onboard soft 



X-ray telescope revealed many jet-like features in which the out-flowing plasma is associated 
with a change in the magnetic field topology of the underlying X-ray emitting structure 
and with the rise of an X-type neutral point location as the reconnection process goes on. 
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Recently, Hinode observations have allowed a more detailed picture of X-jet formation 



and dynamics (Cirtain et al. 2007) and they are suggesting that the heating of the solar 



chromosphere and corona may be related to small-scale ubiquitous reconnection (Shibata 



et al. 2007). Karpen et al. (1998) provided a detailed scenario for the formation of 



current-sheets and their dissipation through reconnection during a chromospheric erupting 
event. They pointed out that magnetic reconnection occurs in a bursty and intermittent 
mode whose unknown physical mechanism is in contrast with usual Sweet-Parker (SP, 
hereafter) and Petschek (PE, hereafter) models of reconnection and it is essential to 
explain the time evolution of flaring events. Their 2. 5- dimensional numerical simulations 
of magnetic reconnection in the lower solar atmosphere showed also the formation of 
two-sided outflows whose dynamics and acceleration should be affected by the vertical 
density variation from the chromosphere to corona's upper layers. These bi-directional jets 
have been observed in the form of simultaneous blue and red Doppler shifts and they are 



characteristic structures present in most of the UV explosive solar events (Aschwanden 



2006, and references therein). 



Large-scale chromosphere/corona models are carried out within MHD. Yet, it is widely 
accepted that in order to obtain a fast reconnection regime in MHD, it is necessary to have a 



localized enhancement of the plasma resistivity (Hautz & Scholer 1987 Scholer 1989 Kliem 



et al. 2000). In fact, it is believed that the description of a plasma as a resistive fluid is not 



sufficient to obtain the fast reconnection rate expected by the models, for instance, of the 
above mentioned structures. On the other hand, the microscopic mechanisms leading to the 
formation of the invoked anomalous resistivity depend on several factors and have not been 



identified conclusively. New insights on this topic have been recently given by Baty et al. 



(2006) who claim to have maintained a PE-like reconnection process in two-dimensional 



time- dependent MHD simulations without the use of a significantly high localized resistivity. 
However, the onset of this process depends still on a non-uniform, even though small. 
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resistivity. Furthermore, Nitta (2007) shows that it is possible to have in MHD a self-similar 



reconnection process with a continuous transition from a slow reconnection regime to a fast 
one by increasing the magnetic Reynolds number. He considers a new series of solutions of 
two-dimensional MHD equations on a wide range of the magnetic Reynolds number, that 
is from about 10 to 2100. The reconnecting structures change according to the regime the 
system undergoes, passing from an X-point configuration, to an X-O-X point, that is two 
X-points with a magnetic islands placed in between. Finally, for the slowest regime at the 
highest Reynolds number the magnetic islands collapse forming a long current-sheet with a 
Y profile. 



On the other hand, Lapenta (2008) observes the spontaneous development of a fast 



reconnection mechanism on macroscopic scales. A two-dimensional current-sheet evolves 
through two different stages, that is a first slow very elongated SP layer is followed by 
a rapid transition to a fast chaotic reconnection process, and the system does not need 
anomalous resistivity nor driven flows to undergo this evolution. This MHD bursty process 
leading to a chaotic flnal state is due to the destabilization of the SP layer via the tearing 



instability as predicted in previous theoretical studies (Furth et al. 1963 Bulanov et al. 



1979 Loureiro et al. 2007). The resulting chaotic reconnection phase is fast and has the 



fundamental property of being independent both from the Lundquist number as well as 
from the dynamic Reynolds number. In order to observe such dynamics, the horizontal size 
of our numerical box must be large with respect to the thickness of the forming SP layer 
and also a high Lundquist number has to be considered. This allows the SP layer to expand 
till the instability driving the system to a chaotic reconnection regime sets in. During this 
phase, an outflow pattern determining a circulation loop feeding-back reconnection regions 
is observed. The net results obtained by the development of such process are that (a) it 
is much faster, developing on scales of the order of the Alfven time, and (b) the areas 
of reconnection become distributed chaotically over a macroscopic region. Furthermore, 
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recent three-dimensional MHD numerical analyses (Kowal et al. 2009) test and confirm 



the investigation on the effects of turbulence on MHD magnetic reconnection (Lazarian & 



Vishniac 


1999; 


Kowal et al. 


2009 



speed of the process is independent both of Ohmic and anomalous resistivity. 

The above described picture does not include the effects of a density gradient either 
parallel or orthogonal to the two-dimensional current-sheet. In compressible conditions, the 
presence of a density variation may influence the onset and the evolution of reconnecting 
field lines and the resulting jet formation and acceleration. By introducing a stratification 
due to a gravitational field, the reconnection process evolves towards a traditional 
two-dimensional pattern passing through an initial stage whose features and duration 



depend on the imposed stratification (Galsgaard & Roussev 2002). In general, a lower 



reconnection rate is found and it decreases according to a faster magnetic field expansion 
which produces a more complex reconnection dynamics. 

The present work has as primary goal to study the spontaneous fast chaotic reconnection 
process in the presence of a strong density gradient modeling the solar atmosphere from 
the high chromosphere up to the low corona. There, the rapid density variation is well 



approximated with a constant plasma pressure (Syrovatskii 1981) where the temperature 
gradient balances the density variation. Hence, the gravitational field can be neglected and 
a force-free configuration of the magnetic field is usually assumed. This reconnection process 
via a chaotic evolution is a promising candidate to explain the cross-scale process linking 
the small-scale reconnection dynamics to the triggering and ejection of large-scale explosive 



phenomena, such as such as solar flares and coronal mass ejections (CME hereafter) (Tajima 



& Shibata 2002). We move from the two-dimensional MHD model of a current-sheet in the 



high chromosphere/ low corona provided by Yokoyama & Shibata (2001) and we determine 
the conditions and the features of a spontaneous transition from a slow reconnection regime 
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to a fast dynamics leading to a chaotic macroscopic state of the reconnection site. 

In the following sections, we introduce the numerical setup, the initial conditions, the 
perturbations, and the parameters for our analysis (Sect. [2]). In Sect [3] and |4| we present 
the results of the simulations for the several cases. We observe the development of a 
two-stage process with the onset of a fast reconnection. The dynamics proceeds with the 
fracture of the initial reconnection site in several structures determining a final turbulent 
state. Moreover, the density jump from the chromosphere to the corona is critical in giving 
an asymmetric development of the global structure, determining an oriented path for the 
overall flow and current- density evolution. We observe also an overall final accelerating jet 
outward to the upper layers of the corona. In Sect. |5} we summarize the results drawing 
our conclusions and pointing to fundamental open questions and future directions. 



We define a longitudinal or stream- wise direction (z), oriented away from the sun, 
wherein open boundary conditions are set at 2; = and L^, and a transverse or cross-stream 
direction (x) along which all the equilibrium quantities vary and where we set reflecting 
boundary conditions. Hence, we solve numerically the viscous-resistive compressible 2.5D 
MHD equations. In dimensionless form we have the following set of equations 



2. Numerical settings 



P 



dp 
di 
dv 

'di 




V ■ pv 




(2) 




A 



u 



(3) 



(4) 
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where p is the mass density, B is the magnetic field, v is the fiuid velocity, X is the specific 
internal energy, J is the current density, p is the fiuid pressure, and 11 is the symmetric 
rate-of-strain tensor defined as 

n = ^ [Vv + Vv^] . (5) 

In Eq. ([IJ-Q, IZx and IZy are the Reynolds numbers measuring respectively the global 
kinematic and dynamic viscosity in the numerical box: in our simulations we always 
consider 1Z\ = IZu and so we refer only to the dynamic Reynolds number IZ^. The quantity 
S is the Lundquist number measuring the global explicit resistivity we set within the 
domain. We use an ideal equation of state with polytropic index 7 = 5/3 such that 
p = p (7 — 1)X. The equations are normalized according to reasonable values for the high 
chromosphere /low solar corona: the number density p/m ^ lO^cm"^; the length scale 
L 1000 — 2000 km, where L = current-sheet width = = Lz/80; the Alfven velocity 

va ~ 400 — 600kms~^ (approximately three times the sound speed). So, we have a time 
scale of tA = L/va ^ 1-5 — 5 s. 



We use the three-dimensional code, FLIP MHD (Brackbill 1991). We consider 600 (in 
x) X 960 (in z) lagrangian markers arrayed initially in a 3 x 3 uniform formation in each of 
the 200 X 320 cells of our numerical grid. This has the following sizes 

xG[-10, 10] withdx = 0.1 

(6) 

;zG[0,80] withd2 = 0.25. 

As shown in Fig. [T| our system consists in a current-sheet determined by a force-free field 
configuration given by the following equations 

Boz{x) = tanhx (7) 
Boy{x) = — sechx, (8) 

which implies the current-sheet width is equal to the normalization length L and it is 
resolved by 30 Lagrangian markers. As already pointed out in the introduction, our 
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computational domain has transversal sizes large enough to allow a complete development 
of the system dynamics. The density is modeled as a step function 



C 



1 



or as a hyperbolic tangent function 

1 + C 

P -- 



+ 



if z < zq = 



if z > zq 



''-^Uanh^"-"^ 



(9) 



2 2 \ a 

where C is the ratio between the density in the region {] < z < L^j'i (modeling the high 
chromosphere) and that one in L^/S < < (the low corona) and we choose it to 



(10) 



range from 10^ up to 10^ as considered in Yokoyama & Shibata (2001). In Eq. (|9|, zi^ 



is the position of the discontinuity resembling the transition region, whereas in Eq. (10) 



it represents the transition region's center; the parameter a allows us to assume different 
thicknesses of the transition region. 

Reconnection is set in by a magnetic field perturbation defined by the out-of-plane 
component of the magnetic vector potential 



Bl = V X {Axyfty 



with 



Axy = e exp [—Kx{x — Xp)'^ — Kz{z — Zp)"^^ cos Kx{x — Xp) cos Kz{z — Zp) (12) 



providing a GEM-like perturbation (Birn et al. 2001) but localized in {x,z) = {xp, Zp) by 



the exponential factor. We referred to the model and simulations of Yokoyama Shibata] 



(2001), but with the intention to study the possibility of a natural-evolving slow-fast 



reconnection transition in the low corona with the development of a final turbulent state. 
So, after having performed a few test simulations, we set the initial perturbation amplitude 
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[e) to 0.5 and the perturbation wave-numbers in the x {kx) and z (k^) direction to tc/L^ 
and Tr/4Lz respectively. In general, we consider that the perturbation is located exactly at 
the center of the x— and axis, that is in the low corona: {xp, Zp) = (0, Lz/2). Other cases 
with different Zp (L^/S, L2/4, 0) but always Xp = have been considered and performed. 



3. Simulation parameters 

3.1. Viscosity and resistivity 

Fig. [2] shows the reconnected flux as a function of the simulation time for several 
different values of the Lundquist and dynamic Reynolds number, both ranging from 10^ 
to 00, even though the actual upper limit is given by the intrinsic numerical effects. The 



plasma (3 is set for these simulation to 0.2 (see § 3.2). As we use a non-selfconsistent 
perturbation that is not an eigenstate of the system, we have a small transient time of 
about a decade of time steps. Afterwards, the system displays the same behavior, except 
for very low S or 71^: it undergoes a two-stage evolution consisting in a slow reconnection 
phase followed by a rapid increase of the reconnected flux which corresponds to a fast 
reconnecting regime. The viscosity and resistivity affect the dynamics of the system in a 
different way. Viscosity is effective only when it is set to high values {TZi, = 10^). Resistivity 
has the main role in determining both the behavior of the system in the first slow phase 
and in switching it to a fast regime. If S is too low, as described by the dotted line in 
the above mentioned figure, the system remains in a steady-state reconnection process 
where the corresponding width of our current-sheet is constant in time until the end of the 
simulation. By increasing S, the slow phase shortens and we observe the triggering of a fast 
regime whose starting point and dynamical features appear to be independent of the value 
of resistivity. If we increase S further until the numerical dissipation dominates {S 00), 
the two-stage evolution continues (as long as viscosity has a negligible value). 
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As long as the reconnection process remains in the first slow reconnection regime, 
the behavior of the system for different values of the resistivity is well-described by the 
SP theory. In fact, the scaling of the reconnection rate, Ai, prescribed by SP theory is 
that Ai = iS^^/^ = r]/{CvA), where 77 is the plasma diffusivity, C the longitudinal size of 
the diffusion region and va is the upstream Alfven speed. In Fig. |3} we show the details 
of the slow regime (as in the small frame of Fig. [2] but) for the curves corresponding to 
S = 10"^ (solid line), S = 10^ (dotted line) and S = 10^ (dashed line). The observed 
reconnection rates are consistent with the expected behavior, even though our basic field 
and perturbation settings does not allow to measure correctly £ and to consider it fixed for 
all the simulations. 

In the subsequent discussion, the focus will be on 5 = 71^ = 10"^. 



3.2. The plasma beta 

As previously pointed out, we consider the simplified assumption, yet reasonable for 
the solar environment, that a temperature gradient is provided in order to balance the 
density variation. Hence, we have an overall uniform and constant initial plasma pressure 

Po = = 2 • (13) 

In Fig. |4], we show the reconnected flux as a function of the simulation time for different 
values of the beta parameter ranging from 0.2 to 2. As the Alfven velocity is fixed, the 
change of (3 implies we are varying the compressibility of our system. In the first slow and 
diffusive SP reconnection regime, the difference in the compressibility does not influence the 
evolution at all. In fact, the dynamics depends on the global system's resistivity/viscosity 
that is fixed by = 5 = lOl Yet, as soon as we observe the triggering of the second 
faster reconnection phase, which corresponds to system's structure evolving on smaller and 
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smaller length-scales (see below), it is evident that for higher compressibility, the slope 
of the reconnected flux becomes steeper and hence the corresponding reconnection rate 
becomes faster. So, the main features of the fast-reconnection phase depend critically on 
the compressive effects in the system and no longer on the global resistivity/viscosity. For 
realistic simulation below, we consider the accepted value (3 = 0.2 to place the current-sheet 



approximately in the high chromosphere/low corona (Yokoyama & Shibata 2001 Gary 
200T|. 



Results and discussion 



The evolution of a current-sheet in the presence of an orthogonal strong density 
variation is characterized by a two-stage evolution whose main features are described in 
Fig. [5| The contour of the out-of-plane component of the current density (Jy) is shown 
superposed on the same component of the magnetic vector potential (solid black lines) at 
four different instants. At about t 18 — 20 (top panel), the system is evolving according 



to the SP theory (see § 3.1) and the slow reconnection regime proceeds with the elongation 
and thinning of the current sheet in the low density region (second panel of Fig. [sj t ~ 36). 



As described in the previous work by Bulanov et al. (1979) and in recently ones (Loureiro 



et al. 2007 Skender & Lapenta 2009), the system reaches a configuration such that the 



ratio between the length and the width of the current sheet makes it tearing-unstable. 
This is evident also in Fig. |6] where diffusion region's aspect ratio is shown as a function 
of the simulation time: at about t ^ 40 a maximum is reached, the instability process 
rapidly evolves and the initially formed SP diffusion region starts to destabilize. Due to 
the disruption of the main current sheet, several structures form each in turn determining 
a reconnection site (third panel from the top of Fig. [sj t ~ 52): those drive the initial 
configuration to smaller and smaller scales (bottom panel). During all the process, the 



- 13 - 



presence of the strong density gradient produces an asymmetric development determining 
an oriented path for the overall flow and current density evolution. As the two-stage 
reconnection proceeds, we observe the slow upward movement of the initial diffusion region 
while it is breaking up into several parts. The islands that form start to move forward and 
backward pushed by the reconnection jets and a chaotic configuration is more and more 
evident both in the most intense current-density region still in fragmentation as well as in 
the macro-structure emitted towards the upper layers of the corona and that resembles the 



turbulent reconnection as proposed by Tajima & Shibata (2002) or Lazarian & Vishniac 



(1999). Furthermore, due to the presence of the density wall the buildup of magnetic flux 
of opposite polarity is taking place and post-disruption arcades are forming in the low 
corona. As the process goes on, converging motions towards the polarity inversion line lead 
to the shrinking of the current-sheet at the base of these growing magnetic structures. This 
situation probably leads to a release of magnetic energy in the solar atmosphere via a new 
reconnection process. 

The onset of a macroscopic circulation pattern linking all the small reconnection sites 
sustains the strong increase of the reconnection rate as seen in Fig. |2] (case with S = 10^). 
The net effect is an overall acceleration of the system that is dragging forward all the 
structures and therefore expelling plasma in the upper corona. In Fig. [7| the contour of the 
out-of-plane magnetic vector potential at t ~ 96 is shown with the velocity field superposed 
and represented by properly scaled arrows: we can observe the well-developed circulation 
pattern during the multiple reconnection process and the effective acceleration reconnection 
jet. The motions driving the post-disruption arcades rising in the low corona are also 
evident. 
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4.1. Reconnection fluxes 

The energy transfers provide information on the structures moving down towards the 
high chromosphere or up towards the upper layers of the solar corona. We consider the flux 
of the Poynting vector, the enthalpy flux and the kinetic flux, respectively defined as 

S = E X B (14) 
H = (T + p)v (15) 
K = ^p|vpv (16) 

through a surface extended from x = — 5 to x = +5 and located at z = 20 (downward flux) 
and z = 80, the higher boundary (upward flux). In Fig. [sj it is possible to observe these 
profiles as a function of time through the two-stage process: the top panel describes the 
upward flux, whereas the bottom panel describes the downward flux. The mass flux through 
the same surfaces is shown as reference. We observe that the fast process determines a 
sensitive rise in the upward mass flux associated with an increase of the kinetic and enthalpy 
fluxes: the multiple chaotically-moving plasmoids accelerate upward and the plasma is 
characterized by an increase of the compressional heating effects. The chaotic fast stage 
shows a initial increase of all fluxes and a weak maximum is reached. Afterwards, a real 
impulsive event occurs and it has a duration of about At ^ 10 that is 15 -i- 50 s according 
to the normalization values provided in § [2] Furthermore, we obtain the following physical 
values for the energy fluxes associated with the impulse 

pv, ^ 4 ^ 0.8 ^2.4 ■ lO^^particle/cm^s (17) 
^ 3 ^ 0.6 ^6.6 ■ 10® erg/cm^ s (18) 
K, ^ 1 ^0.2^2.2- 10® erg/cm^s (19) 



which show a reasonable agreement with observations of solar explosive phenomena (Tsuneta 



1996 Shibata 1996). The Poynting up-flux increases during most of the fast process and it 
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reaches a constant value of about Sz ~ 1.2, that is 0.2 ^ 2.6 • 10® erg/cm^ s, but it drops 
out in the impulsive phase. So, the rapid emission of a highly structured accelerating mass 
towards the upper layers of the corona does not come with a corresponding impulsive release 
of magnetic energy that is instead converted into the other forms of energy. With regards 
to the downward-fluxes, only small variations are observed. This behavior is due to the 
presence of the density gradient that, as already discussed, allows the pile-up of magnetic 
flux with the formation of post-disruption arcades. In Fig. [9} it is possible to observe the 
contour of the specific internal energy at t ~ 96 with respect its initial value, I/I{t = 0). 
It is evident the turbulent configuration of the most dense structures whose internal energy 
rapidly increased during the evolution. In particular, we note the high values obtained on 
the post-disruption arcades in the low corona. 



4.2. Effect of perturbation location 

As already pointed out at the end of § [2} we performed simulations with different values 
of Zp, but no substantial qualitative differences are present as a function of the position of 
the initial perturbation, at least in the cases considered. For instance, if perturbation is 
located on the transition region {xp = and Zp = Lz/8) modeled as a discontinuity between 
the chromosphere and the corona, we observe the same current-sheet dynamics as described 



in the previous sections, even though two clear differences can be observed in Fig. 10, where 
the out-of-plane current density (in color) and magnetic vector potential (solid black lines) 
are shown for this case at t ~ 96: 

• the density step does not allow the current-sheet to develop into the high chromosphere 
at all and this determines that the whole structure is elongated towards upper coronal 
layers. The accelerating jet still presents a chaotic internal structure, but it has a 
wider aperture with respect the previous cases. 
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• the motion of the overall structure towards the upper coronal layer allows again the 
formation of post-disruption arcades that develop slower than they do in the previous 
cases. By the end of the simulation, we are not able to determine any visible sign of a 
starting reconnection process at the base of the corona {Lz ~ 10). 

The upward fluxes in this case do not change with respect the previous ones, while 
the downward fluxes are affected by the presence of the density step determining the 
accumulation of magnetic flux that encompasses the surface through which the flux is 
computed. 

In the case the density is modeled by an hyperbolic tangent, the same dynamics is 
observed as long as a steep profile is considered: the smoother the slope of the density 
function is, the less sensitive the accumulation of magnetic flux, the rise of the main 
reconnection site and the overall acceleration of system is. Finally, with a very low density 



variation, we recover the results obtained by Lapenta (2008), even though the presence of 



open boundary conditions instead of periodic settings alters the detail of the evolution. 



Conclusions 



We presented a 2.5-dimensional model of the current-sheet evolution located in the 
solar atmosphere, in particular across the region where the chromosphere and the corona 
meet. The transition region is modeled as a density and temperature discontinuity or with 
an hyperbolic tangent profile. Previous studies of such systems, most of which focusing on 
triggering events and dynamics of solar flares ( Yokoyama fc Shibata||2001 ), on filaments (Lin 



et al. 2008 ) or on the evolution and observations of current-sheets forming in the aftermath 



of coronal mass ejections (Poletto et al. 2004 Bemporad et al. 2006), always demanded 



the presence of locally enhanced plasma resistivity to explain the rapid increase of the 
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reconnection rate to allow the fast evolution underlying explosive phenomena in the solar 
atmosphere. Yet, in the present paper we demonstrate a natural two-stage evolution of 
a current-sheet with a spontaneous transition from a low to a fast reconnection regime 
without the need to include any anomalous effect and in the presence of a strong density 
gradient. 

We analyzed such mechanism by means of high resolution simulations. We observed 
how the step density jump in the transition region produces an asymmetric development 
of the structure, determining an oriented path for the overall flow and current density 
evolution. Despite the simplifications of our model, the global dynamics reported resembles 
the tearing evolution of solar structures both just before and immediately after observed 
and/or modeled explosive phenomena. Several observational clues can be mentioned to 
support the scenario described in the present paper. For example during solar flares, hard 
X-ray and radio emissions present many properties ascribed to the presence of small-scale. 



fragmented, bursty magnetic reconnections (Aschwanden 2002). Post-CME current-sheets 



are observed to have a much larger thickness than what MHD theory prescribes (Bemporad 



2008). Chaotic magnetic reconnection in a current-sheet could explain this observed 



behavior. 



In the future, an "expanding box" setting (Grappin et al. 1993 Rappazzo et al. 2005) 
of our simulations may be considered in order to follow the evolution of the current-sheet 
in a simpler Cartesian geometry with spherical effects included. This would allow us also 
to properly study the formation and evolution of outward and downward flows as seen 



in Ha observations of highly dynamic filaments and surges in active regions (Lin et al. 



2008). Furthermore, the analysis of the frequency distribution of time scales can possibly 
discriminate between a fractal current-sheet evolution or turbulent reconnection, since 
the former process is scale-free and it generally produces power law distributions, while 



turbulent processes are controlled by incoherent random processes and hence they generally 



give exponential distributions (Aschwanden 2006). Further extension to three-dimensions is 



needed to fully address realistic configurations. 



The research leading to these results has received funding from the European 
Commission's Seventh Framework Programme (FP7/2007-2013) under the grant agreement 
n° 218816 (SOTERIA project, www.soteria-space.eu). The simulations shown were 
conducted using processors on the VIC cluster of the Vlaams Supercomputer Centrum at 
K.U. Leuven (Belgium). 
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Fig. 1. — (color online) Schematic reconstruction of our basic system: it is possible to 
visualize the initial stream-wise component of the magnetic field (in blue and red) and the 
density step at L^/S. 
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Fig. 2. — Big plot: reconnected flux as a function of the simulation time for different values of 
the Lundquist number, S, and the dynamic Reynolds number, TZi,. In all cases but (<S, TZ^) — 
(10^, 10^) (dotted line), we can observe an abrupt transition from a slow reconnection process 
to a fast regime. Small plot: detailed view of the slow SP regime. 
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Fig. 3. — (color online) Reconnected flux as a function of the simulation time but t < 50 (slow 
SP regime) for three cases among those shown in Fig. [2] Reconnected rates are calculated 
and shown. 
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Fig. 4. — Reconnected flux as a function of the simulation time for different values of 
the beta parameter with the Lundquist number and the dynamic Reynolds number set to 
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Fig. 5. — (color online) Contour of the out-of-plane current density, Jy, (in color) superposed 
on the same component of the vector potential (solid black lines) at four different instants 
(from top to bottom): t = 17, 36, 53, 96. Here, the density is modeled as a discontinuity with 
the step located at z = L^/S = 10. For a better visualization, only the central part of the 
numerical box is shown in the images. 
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Fig. 6. — Diffusion region's aspect ratio (length/width) as a function of the simulation time. 
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Fig. 7. — (color online) Contour of the out-of-plane magnetic vector potential at t ~ 96 
(solid colored lines). The velocity field is also shown and it is represented by properly scaled 
arrows. For a better visualization, only the central part of the numerical box is shown in the 
image. 
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Fig. 8. — Top panel: energy and mass fluxes through a surface extended from x = —5 to 
X = +5 and located at the higher boundary, 2; = 80 = (upward fluxes). Bottom panel: 
same as the top panel, but the surface for the computation of the fluxes is located at 2; = 20 
(downward fluxes). 





Fig. 9. — (color online) Contour of the specific internal energy with respect its initial value, 
XjXit = 0), at t = 96. For a better visualization, only the central part of the numerical box 
is shown. 
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Fig. 10. — (color online) Contour of the out-of-plane current density, Jy, (in color) superposed 
on the same component of the vector potential (solid black lines) at t = 96. Here, the 
density is modeled as a discontinuity with the step located aX z = L^/S = 10 and the initial 
perturbation is located at Xp = 0,Zp = Lz/8. For a better visualization, only the central 
part of the numerical box is shown. 



